Environment-assisted quantum transport and trapping in dimers 
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We study the dynamics and trapping of excitations for a dimer with an energy off- set A coupled to an external 
environment. Using a Lindblad quantum master equation approach, we calculate the survival probability II (t) 
of the excitation and define different lifetimes r s of the excitation, corresponding to the duration of the decay 
of n(t) in between two predefined values. We show that it is not possible to always enhance the overall decay 
to the trap. However, it is possible, even for not too small environmental couplings and for values of A of the 
order (9(1), to decrease certain lifetimes r s , leading to faster decay of II (£) in these time intervals: There is an 
optimal environmental coupling, leading to a maximal decay for fixed A. 

PACS numbers: 05.60.Gg, 05.60.Cd, 71.35.-y 



I. INTRODUCTION 

Recent years have seen a growing interest of different com- 
munities in coherent energy transfer. For instance, it was 
noted that in the light-harvesting process of photosynthesis 
quantum mechanical features of the transfer of excitations 
(excitons) created by the incoming solar photons might play 
an important role for the transfer's efficiency 01 0]. Mod- 
els of the (coherent) exciton dynamics in the light-harvesting 
complex of the photo synthetic unit show that the environ- 
ment does not necessarily destroy all coherent features - even 
at ambient temperature - but also can support the coherent 
transfer of excitations [3-7]. In particular, Rebentrost et al. 
have given a detailed analysis of excitonic dynamics in the 
Fenna-Matthews-Olsen complex introducing the concept of 
environment-assisted quantum transport |8]. They showed 
that the maximum transport efficiency is reached for deco- 
herence rates comparable to the difference between the onsite 
and the coupling energies. Cao and Silbey showed that ex- 
citon trapping can be optimized by suitable choices of, e.g., 
decoherence and trapping rates [9fl. 

Rapid experimental advances in ultra-cold gases allow to 
control atoms to a large extent. This offers the possibility 
to study coherent transport and the effect of environmental 
changes (e.g., an increase in the temperature or of an exter- 
nal field). Rydberg gases turn out to be very well suited to 
study the dynamics of excitations II 1 OU 1 211 . At ultra-low tem- 
peratures the dynamics in an ensemble of atoms can be effi- 
ciently modeled by continuous-time quantum walks H 1 311 - By 
properly adjusting specific single Rydberg atoms, the moving 
excitation can be absorbed by these atoms LL5D . An anal- 
ogous process is found in the light-harvesting process, where 
the exciton eventually will reach the reaction center, where the 
exciton 's energy gets absorbed and converted to chemical en- 
ergy. In both cases, the decay of the probability TL(t) not to get 
absorbed monitors the excitation's dynamics. If the process 
is purely coherent, U(t) shows distinct quantum-mechanical 
features such as a power-law decay in certain, experimentally 
relevant time-intervals fl3[ [l4ll. 

Since usually the systems (light-harvesting complexes or 
Rydberg gases) are not isolated from their environment, we 
will model the dynamics by the Lindblad quantum master 
equation (LQME) for the reduced density operator of the sys- 



tem. However, it should be noted that this approach is only 
valid in a limited parameter range of the coupling to the en- 
vironment [16]. Our model system is a dimer, represented by 
two coupled two-level systems, one of which acting as trap. 
We note that various systems, e.g., with radial symmetry and 
a trap in the center [17] or with strong coupling between all 
nodes and weak coupling to the trap, can be mapped onto the 
dimer. 



II. MODEL 

In general, we consider a small system S (network) of N 
nodes coupled to an environment. Each node of the network 
S represents a single two-level system. The basis states \j) 
which are associated with excitations localized at the nodes 
j = 1, . . . , N span the Hilbert space of S alone. The Hamil- 
tonian of the total system, H tot , comprises three parts: the 
Hamiltonians Hs for the network S, Hr for the environment 
(reservoir), and Hsr for the coupling between system and 
environment: if to t = H s + H R + H S r. 

Within a phenomenological approach, the Hamiltonian 
Hs, which incorporates trapping of excitations at the nodes 
m e M, M C {1, . . . , N}, is given by H s = H - iT, 
where Hq is the network Hamiltonian without any trapping 
and iT = iT^ m \m)(m\ (T > 0) is the trapping operator, 
see Ref. 1 14] for details. 

Now, the dynamics of the density operator of the total sys- 
tem is governed by the Liouville-von Neumann equation. In- 
tegrating out the environmental degrees of freedom and as- 
suming bilinear couplings between system and environment 
and the Markov approximation lead to II 1 811 

p(t) = -i[H , p(t)} - {T, p(t)} + V[p{t)\, (1) 

where V[p(t)] is responsible for all decoherence effects. 

Under certain conditons, such as a weak coupling between 
the system and the environment, Eq. (Q]) leads to the Lindblad 
quantum master equation (LQME), where the Lindblad oper- 
ators, acting on the Hilbert space of S, mimick the influence 
of the environment on the dynamics lHHl . Considering only 
localized initial conditions Pfc(O) = \k)(k\ and Lindblad op- 
erators in the form V\\j)(j\ which act only on the diagonal 



2 



elements of p(t), the LQME reads Ut 

p k (t) = -i[H ,p k {t)]-{T,p k (t)} 

N 

-2A^(p,(t)-(j|p,(t)|i))|i)(i|, (2) 

where the initial condition is labeled at the reduced density 
operator by the subscript k. The parameter A represents the 
"strength" of the coupling to the environment. Now, the tran- 
sition probabilities irkj (t) from node j to node k follow from 
the diagonal elements of pk(t), i.e., 7Tkj(t) = (j\pk(t)\j). 
In general, the mean survival probability is defined as 



U(t) = (tr[p k (t)]) h 



1 



N 



N -M 



where (•)& denotes the average over all possible initial nodes, 
i.e., all nodes but the trap nodes. Note, that Eq. © slightly 
deviates from the definition used in |[l4ll for the purely coher- 
ent case; there, only final nodes j ^ m have been considered. 
Equation ® also accounts for the fact that an excitation has a 
probability of going from a node m to a node j ^ m. Thus, 
1 — IL(t) describes the total probability of energy dissipation 
up to time t. 



III. RESULTS 

In order to obtain results which are not blurred by the net- 
works complexity, we consider a dimer model, which is com- 
posed out of two nodes, one of which acts as the trap. How- 
ever, it is also possible to map more complex networks onto 
an effective dimer: One example is a network where the cou- 
plings Vn between the non-trap nodes are very strong and 
only a single node is coupled to the trap with a coupling 
V <C Vn, see Fig. Q] or 111 911. Another example is a ring of 
nodes which are all (weakly) coupled to a trap in the center of 
the ring; if the initial excitation is fully delocalized over the 
ring, the effective dynamics can also be modeled by a dimer: 
the non-trap node of the dimer is given by a superposition of 
all the nodes in the ring and the trap of the dimer is identified 
with the trap in the center of the ring. 
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FIG. 1: Schematic illustration of the dimer model. 

Now, consider a dimer coupled to an external bath. The 
Hamiltonian of the dimer without any coupling to the sur- 
roundings can be expressed in matrix notation by 



Hq = Ho 



iT 



E 1 - 
-V E 2 



V 



iT 



(4) 



where E\ and E 2 are the on-site energies and V is the cou- 
pling between the two nodes. Note that in the following 
we will express all varying parameters in units of V. The 
eigenvalues of H are E± = E\ ± Ve ±( ^ A , where 0a = 
arcsinh[(A - iT)/2V] and A = \Ei - E 2 \, and where we 
assumed T < 2 V, such that (for A = 0) the dimer is not over- 
damped 111 611 . The bi-orthonormalized eigenstates of H are of 
the form 



i 



V2cosh</. A V±e T0A/2 



5 ±<M/2 



(5) 



and 



^cosh^ V ±e±0l/2 



(6) 



where (j)\ is the complex conjugate of </>a- In the limit A — » 
we recover the results of 11611 . Note also, that our model is 
similar to the one studied by Cao and Silbey in |2D. There, the 
authors obtained analytical estimates for the (quantum) mean 
first passage time. 




FIG. 2: 

for A = 
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(Color online) Survival probability II (t) for T — V — 1: (a) 
= and different A, see Eq. (0. (b) solution of the LQME for 
1 and different A. (c) same as (b) but for A = 1.5. 
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A. Survival probabilities 



t\. decrease of II (t) from 1 to e 1 , 



Let the initial excitation now be placed on the non-trap node 
1, i.e., (l|pi(0)|l) = 1. FromEq. ©, the survival probability 
follows as U(t) = (l|pi(t)|l) + (2|pi(t)|2). Without de- 
phasing (A = 0) the survival probability is obtained by direct 
diagonalization of Hs, yielding the analytical result 



-Ft 



n(t) 



|cosh(<^ 



|cosh(*-/(r,A,V) + <A A )r 



|sinh(t-/(r,A,y))r 



(7) 



with/(r,A,y) = 

r = y/y 2 + 4r 2 A 2 and y 
A -> yields d]: 



y)/2 + iy/(r + y)/2j, where 
: 4V 2 - T 2 + A 2 . The limit 



H(t) 



-rt 



cos^ 



cos 2 (Vt cos(0 r )— (/) r )+sin 2 (yt cos(</> r )) 



(8) 

with (j) r = arcsin(r/2F). Letting now also r —> 0, one ar- 
rives at the usual resonant Rabi oscillations for the dimer Il20ll. 

Figure[2a) shows H(t) for A = and for different values of 
A: Increasing A leads to a slower decay of II (t). While small 
values A = 0.1 do not lead to strong deviations from the 
A = results, larger values of A clearly shift the whole H(t) 
curve upwards. This can be interpreted as the rudimentary 
onset of localization, the disorder being the off-set A. Note 
also that all curves have a plateau-like region where the slope 
decreases. By increasing A this "plateau" is shifted to larger 
values of H(t). This plateau is due to the fact that the dimer 
only has two nodes. For larger networks, the averaging in 
II(t) diminishes such effects. 

Now, consider the case with a small trapping strength (T <C 
V) and let also the coupling be small (A « 7). By combining 
the results for A = and for T = , and by expanding all 
terms except exponentials to first order in T and A, results in 
the simple exponential decay U(t) ~ e~ rt for not too short 
times. 

Plots of H(t) for different values of A ^ and A ^ are 
shown in Fig. [2b) and (c). The behavior is similar both for 
small values A = 0.1 and larger values A = 1.5: An increase 
of the coupling A leads always to a slower decay of H(t). 
Moreover, for fixed A, increasing A also leads to a slower 
decay of II(t), compare, e.g., the solid curves in Fig. Ob) and 
(c). However, both for small and for larger A, the curves do 
not stay equidistant at all times. Therefore, on different time 
scales the trapping efficiency, here defined as the time it takes 
to decrease IL(t) from one value to another, can vary. 



B. Lifetimes 

In order to obtain a quantitative measure for the efficiency 
of the transport from node 1 to the trap node 2, we define dif- 
ferent lifetimes of the excitation, based on the decay of U(t) 
between two given values: 



T2'. decrease of H(t) from e _1 to e -2 , 

T3: decrease of U(t) from e -2 to e -3 . 

see also Fig. [3] Thus, after t\ = t\ about 63% of the proba- 
bility has been absorbed by the trap, after t 2 = Tl + r 2 about 
86%, and after £3 = n + r 2 + T3 about 95%. By consider- 
ing the lifetimes r s individually, we are able to quantify the 
rate l/r s at which a certain amount of probability gets to be 
absorbed given the initial value of U(t s -i). 
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FIG. 3: Definition of the lifetimes n, T2, and T3, see text. 

Additionally, we also consider the time r^, which is ob- 
tained by rewriting the LQME, Eq. ©, in the form pk(t) = 
Cpk(t), whose formal solution is Pk(t) = exp(Ct)pk(0), 
where p k (t) = ((l|p fc (*)|l), (l|p*(t)|2), . . . , (N\p k (t)\N)) 
is an TV 2 -dimensional vector. Diagonalization of C then leads 
to pk(t) = Qexp(At)Q _1 pfe(0), where Q is the eigenvector 
matrix of C and A a diagonal matrix whose elements are the 
eigenvalues of C. The long-time decay rate 1/tqo is then given 
by the smallest of all the real parts of the eigenvalues of C. We 
checked numerically that in our case J °° dt(n\pk (t)\n) is 
identical to the results of Cao and Silbey for the mean first 
passage time, see Eq. (4) in [9] . We note, however, that here 
we focus on the different lifetimes defined by the mean sur- 
vival probability and not on averaged quantities as the mean 
first passage time. 

Figures Hta)-(d) show contour plots of n, T2, T3, and r^, 
respectively, as functions of A and A for V = T = 1. All 
figures show a similar behavior: 

(i) For fixed A and increasing A all lifetimes increase 
monotonically (except rs for small A, see below). 

(ii) For not too small values of A, however, all lifetimes 
except r\ (see below) first decrease with increasing A to 
a minimum and then increase again. 

The latter behavior is exemplified in Fig.Hfe) for A = 1.5. 
Clearly there is a minimum for T2, T3, and at about the 
same value of A « 0.5. Therefore, there is an optimal value 
of A =^ which leads to larger rates l/r 3 (s ^ 1) compared to 
the ones for A = 0. This implies, for instance, that for T2 an 
additional 23% of the probability gets to be absorbed by the 
trap faster if the dimer is coupled to the environment. How- 
ever, it does not neccessarily imply that 86% of the probability 
transferred to the trap after time t 2 gets to be absorbed faster. 
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FIG. 4: (Color online) (a)-(d) Contour plots of the lifetimes n, T2, 
T3, and Too, respectively, as functions of A and A for V — T — 1. 
The arrows at A = 1.5 show the location of the cuts for fixed A and 
increasing A displayed panel (e). 



spread the excitation evenly over the dimer. However, with in- 
creasing values of A one approaches the quantum Zeno limit 
which also leads to a slower decay of the survival probability 
and thus to larger lifetimes. Therefore, the minimal lifetimes 
can be viewed as being caused by an optimal decoherence rate 
A destroying localization due to A. We note that our results 
are in agreement with the findings of Cao and Silbey | 9] and 
with the ones of Rebentrost et al. |8fl. 

The fact that the lifetimes r s (s > 2) lie on the same curve 
for A > 0.5 translates to an (on average) exponential decrease 
of the survival probability: Let the lifetimes r s be obtained 
piecewise, approximately via H(t s -i) — H(t s ) ~ exp(— ar s ), 
where a is some constant. Now, if r s = r for all s, we ob- 
tain U(t s - r ) — H(t s ) ~ exp(— arr), with r = 1, . . . , s — 1. 
Allowing also r = s, one has U(t s ) ~ II (to) — exp(— at 3 ), 
where t s = sr. 

Finally, we remark on one peculiar feature of the behavior 
of T3. For small values of A < 0.1, the lifetime rs has a 
local maximum at A « 0.7, visible as a small "lagoon" in 
Fig-fflc). This is due to the plateau of II(t), mentioned above. 
By increasing A for fixed (small) values of A the plateau is 
shifted to larger values of U(t) and, therefore, can eventually 
lie in the region between II (t) = e _1 and II(t) = e -2 . This 
then might lead to longer lifetimes rs . The other lifetimes T2 
and are not affected by this. 



IV. CONCLUSIONS 



The same reasoning holds also for rs and . The initial ab- 
sorption at the rate 1/n cannot be enhanced for these values 
of A of the order (9(1) by the coupling to the environment. 
This would mean that, in the interval t\ , (the later part of) a 
U(t) curve for A' > A would lie below the curve for A, which 
within our LQME approach is not the case, see Figs.|2](b) and 
(c). However, increasing A to larger values also leads to a 
shallow minimum for r\ as a function of A, see Fig. Ufa) for 
A « 2. 

The increase of the rates also depends on the value of A. 
While an increase of r^ is visible even for small A w 0.1, 
this is not the case for T2 and T3. Thus, in order to obtain a 
sufficiently large increase of the rates T2 and T3, one needs 
to have a value of A of the order (9(1). Intuitively, one ob- 
serves a "competition" between localization and decoherence 
effects: While the energy offset A tends to bind the excita- 
tion at the initial node, the environmental coupling A tends to 



In conclusion, we have studied the trapping of excitations 
for a dimer with an energy off- set and with environmental 
coupling. The survival probability, II (t), not to get trapped 
(obtained from a Lindblad quantum master equation) shows 
distinct features depending on the strength A of the coupling 
to the environment and on the energy off-set A: While it is 
not possible to always enhance the overall efficiency of the 
decay to the trap, it is possible to increase the rate of trapping 
in certain intervals of time, leading to faster decay of proba- 
bility in these intervals. The fastest decay is obtained for an 
optimal value of A which is about half the value of the cou- 
pling V between the nodes. However, a substantial increase is 
only obtained for values of A of the order (9(1). 
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